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Abstract. Space-time evolution and subsequent particle production from minimally viscous 
(77/5=0.08) QGP fluid is studied using the 2nd order Israel-Stewart's theory of dissipative 
relativistic fluid. Compared to ideal fluid, energy density or temperature evolves slowly in 
viscous dynamics. Particle yield at high pr is increased. Elliptic flow on the other hand 
decreases in viscous dynamics. Minimally viscous QGP fluid found to be consistent with a large 
number of experimental data. 

I) . 

;-C ■ 1. Introduction 

Success of ideal hydrodynamics [1], in explaining bulk of the data in Au-|-Au collisions at RHIC, 
CN ■ has led to a paradigm that in Au-|-Au collisions, a nearly perfect fluid is created. However, 

the paradigm of perfect fluid, produced in Au-|-Au collisions at RHIC, need to be clarified. 
{^ \ Experimental data do show deviations from ideal behavior; at large px > 1-5 GeV, in peripheral 

Tij" ■ collisions or at forward rapidity [2j , presumably due to increasingly important role of dissipative 

effects. ADS/CFT correspondence [3] also suggests that in a strongly coupled matter, the 
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'nJ" ' shear viscosity to entropy ratio is bounded from the lower side, ry/s < 1/4tt. At the minimum, 

rj/s = l/47r ~ 0.08 for QGP fluid. Hydrodynamics is applicable only when, rj/s « Tt [4J. In 

^^ \ Au-|-Au collisions at RHIC, at the early time, Tj=0. 6-1.0 fm, temperature is T, ~300-350 MeV, 

which limit the viscosity to rj/ s « 0.9-1.8. ADS/CFT bound is smaller by factor of 10-20 from 
the limiting viscosity. 

K> I Though the theories of dissipative hydrodynamics O El [7] has been known for more than 

j^ ■ 30 years , significant progress toward its numerical implementation has only been made very 

recently [H El El [lOl [HI [I2l [I3l [HI [151 [K] ■ Several groups have solved the equations for causal 
hydrodynamics. The program is still ongoing and consensus is not reached between different 
groups. In the following, we discus some aspect of viscous dynamics and associated results. 

2. Causal dissipative hydrodynamics 

We consider QGP fluid in the central rapidity region, with zero net baryon density and chemical 
potential, ub = 0, fiB = 0. We neglect the effects of heat conduction (fj^B = 0) and bulk viscosity 
and account only for the shear viscosity. We work in the Landau-Lifshitz energy frame. The 
energy-momentum tensor, including the shear stress tensor vr'^'^, is written as, 

T^-^ = {£ + p)ui'u'' - pg^"" + TT^", (1) 



where e is the energy density, p is the hydrostatic pressure, and u is the hydrodynamic 4- velocity, 
normahzed as u^u^ = 1. T^^ satisfies the energy- momentum conservation law. 

In the Israel-Stewart's 2nd order theory of dissipative fluids [7], the dissipative fluxes are 
treated as extended thermodynamic variables which satisfy relaxation equations. The relaxation 
equation for the shear stress tensor read 

Z)7r^'^ = — (7r^'^-2r?V<^u^>), (3) 

where D = u^d^ is the convective time derivative, r] is the shear viscosity coefficient and r^ is 
the relaxation time. V^^u^' is a symmetric, traceless tensor, 

V^i^u"'^ = \[V^'u'' + V^n^] - \A^"'d^u'', (4) 

where V^ = d^ — u^D is the transverse gradient operator and A^'^ = g^^^ — u^u'^ is the projector 
orthogonal to the flow velocity u^. In the time scale r^r, viscous pressure relaxes to 1st order 
value, 2r]'V^^u'^' E]- The viscous pressure tensor vr^*^ is symmetric {ir'^'^ = tt^^), traceless 
(■^u — 0) and transverse to the hydrodynamic 4-velocity, u^ir^'^ = 0. It has 5 independent 
components. There are 10 unknowns (e, p, three components of the hydrodynamic velocity 
u, and 5 viscous pressure components) and 9 equations (4 energy-momentum conservation 
equations and 5 transport equations for the independent components of tt^'^). The set of 
equations is closed by the equation of state p = p{e). We note that presently there is 
disagreement about the form of the relaxation equation as given in Eql3l Some authors |141 [T5] 
have solved the relaxation equation containing an extra term R = [tt'^Tr'^^ + u^7r'^^]Dux. The 
term ensures that throughout the evolution shear stress tensor remains traceless and transverse 
to fluid velocity. Israel-Stewart developed the theory on gradient expansion of entropy density. 
Gradients of equilibrium thermodynamical variables are assumed to be small. Since dissipative 
flows are small, terms like n^^Duy is neglected. Consequently, Eql3] is restricted to situations 
where gradient of velocity is small. In the present paper, we have limited our study to minimally 
viscous fluid (ry/s=0.08). It will be shown later that for minimally viscous fluid, contribution of 
the term R is small and fluid evolution is hardly changed whether the term R is included or not 
in the relaxation equation. 

Assuming longitudinal boost-invariance, we solve Eqsl2] and [3] in (r = \Jt'^ — z'^,x,y,r] = 
2 In j^) coordinates. The assumption reduces the number of unknowns to 6 (e, two components 
of u and three components of shear stress tensor. We choose vr^^, vr^^ and vr^^ as the independent 
shear stress tensor components. There can be other choices also |14tll5j. Six partial differential 
equations are simultaneously solved with the code " AZHYDRO-KOLKATA" developed at the 
Cyclotron Centre, Kolkata. We use the standard initial conditions and equation of state EOS-Q 
(incorporating 1st order phase transition at rc=164 MeV), described in [l]. At initial time 
Ti=0.6 fm, the QGP fluid was initialized with central entropy density Sini = 110fm~^, with a 
Glauber model transverse density proflle. It corresponds to peak energy density ~ 35GeV/ fm^. 
Initial fluid velocity is zero. In viscous hydrodynamics, additionally, one has to initialize the 
shear stress-tensor. Viscous effects are enhanced if initially n^"^ is non-zero rather than zero. 
We choose to maximise the effect of viscosity. Instead of choosing any arbitrary non-zero value, 
we initialise shear stress-tensor to boost-invariant value, vr^^ = vr™ = 2r//3Tj, vr^^ = 0. For the 
relaxation time, we use the kinetic theory approximation for a Boltzmann gas, Tn = j^ ^ -^- For 
the viscosity, we have used the ADS/CFT lower bound r]/s = 0.08. Throughout the evolution, 
bound on the shear viscosity to entropy ratio is maintained. 
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Figure 1. (color on line). For testing the computer code AZHYDRO-KOLKATA, we have 
considered fluid evolution in QGP phase only. In the left panel, contour plot of energy density 
show that the evolution is stable with respect to change in integration step- lengths, maintain any 
given symmetry and do not lead to unphysical maxima or minima. The middle panel shows that 
the fluid at the centre is least affected by the transverse expansion, and closely follow the 0+1 
dimensional causal flow. In both left and middle panel viscosity to entropy ratio is r//s=0.08. 
The right panel shows that ideal hydrodynamic results are recovered when viscosity gradually 
reduces to zero. 
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Figure 2. Evolution of energy density 
in ideal and minimally viscous flow in 
b=0 Au+Au collision. 



Figure 3. Evolution of stress tensor 
TT^'^' and iryy in b=0 Au+Au collision. 
?7/s=0.08. 



AZHYDRO-KOLKATA was tested extensively [13] for numerical accuracy and consistency 
(see Fig[l|. Many details of fluid evolution in minimally viscous hydrodynamics can be found 
in |13) . For completeness, in Fig 121 we have compared evolution of minimally viscous fluid 
and ideal fluid, in a b=0 Au+Au collision, both initialised similarly. As expected, fluid cools 
slower with viscosity. In FiglSl evolution of shear stress tensor components, vr^^ and tt^^ in 
x-y plane is shown. Initially vr^^ and vr^^ both have similar distribution. After a few fm of 
evolution, the distribution starts to differ, tt^^ and vr^^ are related by x ^ y and y — > rr. The 
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Figure 4. (color online) Evolution of 
energy density in a b=0 fm Au+Au colli- 
sion. The black lines are constant energy 
density contours in r(= \/ x^ + y^) — 
r plane in ideal dynamics. The red 
lines are for minimally (r//s=0.08) vis- 
cous fluid with Israel-Stewart's relax- 
ation equation EqHl which neglect the 
term R = [u^'TT''^ + u''7r''^]Dux. The blue 
lines are energy density contours when 
the relaxation equation contain the term 
R. At early time, fluid evolution, with 
or without the term R remain essentially 
unchanged. Only at late time, evolution 
is marginally changed if the term is in- 
cluded in the relaxation equation. 



relation is maintained throughout the evolution. In FiglH black and red lines show the constant 
energy density contours in r — r plane for ideal and viscous fluid respectively. In [141 [TS] it 
was observed that at late time, centre of the fluid cools faster in viscous evolution than in ideal 
evolution. However, we do find a contrary result, even at late time, centre of the fluid cools 
slower than ideal fluid. However, in peripheral region and at late times, we do find that viscous 
fluid cools faster than ideal fluid. As mentioned earlier, relaxation equation Eql3] neglect the 
term R = [n^vr'^ +u'^it'^ ]Dux. To check that indeed the term R contributes negligibly, in FigH] 
we have shown results obtained with the term R included in the relaxation equation. The blue 
lines in FigS] indicate that fluid evolution with or without the term R is nearly identical. For 



minimally viscous fluid neglect of the term R 
is justified. 



u^^ir 
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+ u'^tt'^ ]Du\ in the relaxation equation 



3. Particle spectra and elliptic flow 

Hydrodynamic equations give the space-time evolution of thermodynamical quantities e.g. 
energy density (or temperature), fluid velocity and the shear stress tensor. The information 
is converted into particle spectra using the standard Cooper-Frye prescription. In Cooper-Frye 
prescription, invariant distribution of a particle, say tt", is obtained as. 



dN 



dS^p^/(o)(x,p)[l + </)(x,p)], 



(5) 



dydPpT JT. 

where S^ is the freeze-out hyper surface and f^^'{x,p) is the equilibrium one-body distribution 
function and 0(x,p) = 2T'^LI ) « li is the non-equilibrium correction to the equilibrium 
distribution function. Accordingly, in viscous dynamics, invariant distribution has two parts. 
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dN'"! 



dyd'^pj 
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^j^neq 



dyd^PT dyd'^pT ' 



(6) 



where ^..^^ is the non-equilibrium correction to the equilibrium distribution -r—p, — . Similarly, 



in viscous hydrodynamics, elliptic flow of can also be decomposed in to two parts [Ij, 
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Figure 5. Ratio of non- Figure 6. px spectra of vr" Figure 7. Elliptic flow in 
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are the equilibrium flow (disregarding (p{x,p)) and its non-equilibrium 



where ^2'' and 
correction respectively. 

Since the non-equilibrium correction to equilibrium distribution (p{x,p) « 1 it is then 
necessary that the ratio R = -^-^ — /tIF- — ^"^ -'■■ ^^ Figl5l pT dependence of the ratio R in 
a b=6.5 fm Au-|-Au collision is shown. R depend sensitively on the freeze-out temperature, 
decreasing with lowering Tp. pT range over which hydrodynamics remain applicable also 
increases as Tp is lowered. 

In Figl6] and \7\ we have demonstrated the effect of viscosity on pT spectra and elliptic flow. 
For both ideal and minimally viscous fluid, vr" yield in a b=7 fm Au-|-Au collision is calculated 
from freeze-out surface at temperature Tp=130 MeV. Compared to ideal dynamics particle yield 
is increased in viscous dynamics, more at large px, e.g. at pr=3 GeV, 7r~ yield is increased by a 
factor of ~10. Elliptic flow on the other hand decreases in viscous dynamics. As shown in FiglTl 
while in ideal hydrodynamics elliptic flow continue to increase with pT, it nearly saturate in 
viscous hydrodynamics. It is well known that ideal hydrodynamics under-predict the pT spectra 
and over-predict the elliptic flow. Viscous dynamics seems to remedy the drawbacks of ideal 
dynamics. 

Data fitting is a complex process in viscous dynamics. Apart from initial time, initial energy 
density and fluid velocity (as in ideal hydrodynamics), in a minimally viscous fluid, particle 
production depends on the initial shear stress tensor, the relaxation time (see FigJH] and E]) . 
Viscous effects are enhanced with increasing relaxation time. Effects are also increased with 
non-zero initial tt^*^, a result contrary to [K] where pT spectra or V2 show little sensitivity to 
the initial shear stress tensor. 

Presently we fit the PHENIX data [T7] on differential elliptic flow in 16-23% centrality Au-|-Au 
collisions, by varying only Tp. The other parameters are kept fixed, rj=0.6 fm, Smi = IWfm^'^, 
Vx = Vy=Q, ir^^ = T^yy = 2rj/3Ti,Tr^y = 0, T,r = (irj/sT. In FigfTOl elliptic flow in a b=6.5 fm 
Au-|-Au collision is compared with the PHENIX data for a range of freeze-out temperatures. 
We have shown the equilibrium flow Wg"^ (the dash-dotted lines), the non-equilibrium correction 
^2*^'^ (the dash-dot-dotted lines) and total flow V2 (the solid lines) separately. It is interesting to 
note that ^2'' change marginally from Tp=130-160 MeV . It indicates that equilibrium flow is 
early time phenomena. Most of the flow is generated by the time fluid cools to T=160 MeV, and 
further evolution do not generate significant flow. Non-equilibrium correction is negative and 
as expected, reduces with lowering Tp. The black arrows in FigdO] indicate the pp above which 
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Figure 8. (a) pT spectra and (b) elliptic 
flow in b=6.5 fm Au+Au collisions. Solid 
and dashed lines are for initial 'n^^ = and 
vr^'^ = 2r]/'iTi respectively. 




Figure 10. PHENIX data on the 
Pt dependence of elliptic flow in 16- 
23% Au+Au centrality collisions are 
compared with minimally viscous fluid 
evolution (see text). 
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Figure 9. (a) pT spectra and (b) elliptic 
flow in b=6.5 fm Au+Au collisions. The 
solid, dashed and short dashed liens are for 
TTr=3r]/sT, Qrj/sT and 9r]/sT respectively. 
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Figure 11. Predictions from minimally 
viscous hydrodynamics for V2 in different 
centrality ranges of collisions are com- 
pared with the STAR and PHENIX data 
(see text). 



viscous hydrodynamics break down. As seen in FigllOl for Tp=130 MeV, the total elliptic flow 
agree well with the PHENIX experiment. For comparison, in FigllOlwe have also shown the V2 
in ideal dynamics (the blue lines). In ideal dynamics, comparable fit is not obtained. 

As shown in Figllll minimallv viscous hydrodynamics, with Tp=130 MeV, also reasonably 
well explain the PHENIX data [H] on differential V2 in 0-10%, 10-20%, 20-30% and 30-40% 
centrality Au+Au collisions (the colored symbols) and the STAR data [19] (open symbols) 
for minimum bias V2- The black lines (from bottom to top) are V2 from minimally viscous 
hydrodynamics in b=3.2, 5.7, 7.4 and 8.7 fm Au+Au collisions. They roughly corresponds 
to 0-10%, 10-20%, 20-30% and 30-40% centrality Au+Au collisions. V2 in 10-20% or in 20- 
30% centrality collisions are reasonably well explained, but V2 is under predicted in 0-10% 
centrality collision and over-predicted in 30-40% centrality collisions. Interestingly, STAR data 
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Figure 12. (color online) PHENIX data [20j on vr", K~^ and proton pT spectra, in different 
centrality ranges of Au+Au collisions. Solid lines are predictions from minimally viscous 
hydrodynamics . 



on minimum bias V2 is correctly explained (the blue line in Fig llip . In minimum bias, all the 
centrality ranges of collisions are included. Two opposing effects {v2 under-predicted in central 
collisions and over predicted in peripheral collisions) are cancelled in minimum bias V2- 

Minimally viscous hydrodynamics with freeze-out temperature Tp=130 MeV also reproduces 
the centrality dependence of pr spectra of identified particles . In Fig{T2] we have compared 
viscous hydrodynamics predictions with PHENIX data [20] on vr", K^ and proton p^^-spectra 
in 0-5%, 5-10%, 10-20%, 20-30%, 30-40% and 40-50% centrality Au+Au colhsions. The spectra 
are are normalised by a factor A^ = 1.4, which is reasonable considering that we have neglected 
resonance contribution. While ideal hydrodynamics can reproduce pT spectra only up to 
Pt ^ 1-5 GeV p!], viscous hydrodynamics can reproduce the spectra throughout the pT range. 
Viscous hydrodynamics also reproduces centrality dependence of charged particle multiplicity, 
mean px and pr integrated V2 (see Fig JT^ . 

The present study indicate that minimally viscous QGP fluid, with central energy density 
~ 35GeV/fm^, thermalised at rj=0.6 is consistent with experiments if the hadronic fluid freeze- 
out at ri?=130 MeV. However, it must be mentioned that initial or final condition of the fluid 
may be changed if a different r^r or initial n'^'^ is used. A comprehensive analysis, exploring all 
the variables upon which the fluid evolution depend is required to identify the initial condition 
of the fluid produced in Au+Au collisions at RHIC. Apparently, present simulation contradicts 
simulations in |141 115|. Different initial conditions, choice of system, evolution parameter etc. 
makes direct comparison between different simulations difficult. Song and Heinz [15] studied 
viscous flow in Cu+Cu collisions. They did not compare with data but V2 in Cu+Cu indicate 
that possibly minimal viscosity will be inconsistent with the RHIC data. Romatschke et al |14) 
studied Au+Au collisions. Minimum bias elliptic flow is under-predicted for minimal viscosity, 
data require r]/s ~0.03. However, they used freeze-out temperature Tp=150 MeV and we expect 
they will find better agreement with data if the freeze-out temperature is lowered to Tp=130 
MeV. 



4. Summary 

To summarise, we have solved Israel-Stewart's 2nd order theory to study evolution of minimally 
viscous QGP fluid and subsequent particle production. In viscous dynamics, energy density 
or temperature of the fluid evolve slowly than in ideal fluid, for a fixed freeze-out temperature 
lifetime of the fluid is increased. Particle production is enhanced in viscous dynamics, more at 
large pT- The elliptic flow on the other hand is reduced. Within certain approximation (non- 
zero initial shear stress tensor, kinetic theory approximation for relaxation time etc.), ADS/CFT 
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Figure 13. (color online) PHENIX 
data 120] on the centrality dependence 
of dNch/dy, <PT> and PHOBOS data 
[2T] on pT^-integrated V2 are compared 
with predictions from minimally viscous 
dynamics (the solid lines) and ideal dy- 
namics (the dashed lines). Freeze-out 
temperature is Tp=100 MeV for ideal 
and Tp=130 MeV for viscous fluid. Vis- 
cous hydrodynamics give nearly equiva- 
lent description of centrality dependence 
of dNch/dy as in ideal dynamics. For 
Npart >100, centrality dependence of 
< Pt > is better explained in viscous 
dynamics than in ideal hydrodynamics. 
PHOBOS data on integrated elliptic flow 
is also well explained. 



lower bound on viscosity appears to be consistent with a large number RHIC data in Au-|-Au 
collisions. 
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